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Collective motion is one of the most ubiquitous behaviours displayed by social organisms and has 
led to the development of numerous models. Recent advances in the understanding of sensory system 
and information processing by animals impel to revise classical assumptions made in decisional 
algorithms. In this context, we present a new model describing the three dimensional visual sensory 
system of hsh that adjust their trajectory according to their perception held. Furthermore, we 
introduce a new stochastic process based on a probability distribution function to move in targeted 
directions rather than on a summation of inhuential vectors as it is classically assumed by most 
models. We show that this model can spontaneously transits from consensus to choice. In parallel, 
we present experimental results of zebrahsh (alone or in group of 10) swimming in both homogeneous 
and heterogeneous environments. We use these experimental data to set the parameter values of our 
model and show that this perception-based approach can simulate the collective motion of species 
showing cohesive behaviour in heterogeneous environments. Finally, we discuss the advances of this 
multilayer model and its possible outcomes in biological, physical and robotic sciences. 
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I. INTRODUCTION 
A. Modeling collective motion 

Collective motion is one of the most ubiquitous collec¬ 
tive behaviour displayed by interacting organisms such 
as cells m, bacteria HHZ], invertebrates (in locusts: 
HHini ; in ants : [in US; in honeybees : m) and ver¬ 
tebrates species (in birds: pAlUb] : in fish: p!7lll9] : in 
mammals: [20]) including humans [2TJ[22|. A growing in¬ 
terest to decipher the link between individual behaviours 
and collective patterns has arisen out of these numerous 
observations and led to the development of models sim¬ 
ulating agents performing collective movement inspired 
from birds, mammals or fish, the latter being the focus 
of this paper. 

Different types of model exist to describe fish school¬ 
ing. In self-propelled particles (SPP) models (syn¬ 
chronous [23l [24| or asynchronous Haile]) firstly used 
for computer animation m, interactions between fish 
are mostly composed by a collision avoidance compo¬ 
nent, a alignement component and a cohesion compo¬ 
nent [28j|29]. Similarly, in social forces (SF) models fish 
are considered as Newtonian particles subjected to social 
and physical forces that respectively ensure the cohesion 
of the group and reflect the interaction (drag for exam¬ 
ple) with the environment [30l [31]. Both SPP and SF 
models have inspired several studies in statistical physics 
that aim to characterize features of a large number of 
interacting agents at the collective level [32H36] . Finally, 
kinematic models describe the evolution of the trajec- 
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tory of fish by a stochastic differential equation [37H4Q] . 
This modelling approach has successfully described the 
movement of fish in closed environment and is a con¬ 
tinuous time formulation of a particular case of random 
walks (RW). Random walks describe stochastic trajecto¬ 
ries build by successive random steps that can be drawn 
from a uniform distribution (unbiased random walk) or a 
non-uniform distribution (biased random walked). These 
probabilistic models have a wide range of application, 
from simulating the brownian motion of particles to the 
exploratory patterns of many species including humans 
1311 [42]. 


B. Information perception 

In all mentioned models of schooling, the agents move 
according to their conspecifics (position, speed, orienta¬ 
tion or a combinaison). Multiple hypothesis have been 
proposed to calculate the subset of individuals that in¬ 
fluence the motion of a focal fish: the metric percep¬ 
tion includes all individuals situated within a defined dis¬ 
tance; the topological perception includes the proxi¬ 
mal neighbors; the perception based on Voronoi tessella¬ 
tion includes the fish connected to the focal fish by the 
Delaunay triangulation. However, while these hypothesis 
produce coherent movement of simulated agents, they are 
not sufficiently constrained by known biology. Therefore, 
recent works are now based on visual perception [43ll45] . 
In these models, the focal fish does not interact with its 
neighbors according to their Cartesian coordinates but 
according to their representation in its visual field. 

Thus, theoretical and experimental studies highlighted 
that the visual sensory system of fish is determinant 
for information transfer during collective motion. In¬ 
deed, the comparison of interaction mechanisms showed 
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that a vision-based model outperforms others mecha¬ 
nisms (metric, topologic, Voronoi) in explaining exper¬ 
imental data [45]. In parallel, the increasing knowledge 
on the visual system of fish allow to develop model based 
on sensory systems more coherent with the biological re¬ 
ality. The characterization of the vision of cyprinids re¬ 
veals that they have a wide visual coverage in the hor¬ 
izontal and vertical planes and an acute vision in the 
front-dorsal region [46]. To make a step towards more 
realistic sensory systems, we introduce in this paper a 
new mechanism based on a the perception of 3D stimuli 
in the visual field. We simulate fish-agents that perceive 
stimuli (congeners, environment) according to the solid 
angle -an analogous of the planar angle but in 3D- that 
they intercept in their visual field. 


Metric Topological Visual 
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C. Information processing 

Once the focal fish has perceived potential stimuli, this 
information has to be processed and translated into a 
movement of the individual. In SPP and SF models con¬ 
sidering only the influence of congeners, a vector of inter¬ 
action is computed with all neighbors situated in a delim¬ 
ited range (metric models) or according to their proxim¬ 
ity rank (topological models). Then, the focal fish moves 
along the resulting force computed by a weighted sum¬ 
mation process of the different interaction forces applied 
on the focal fish. While it results in coordinated motion 
of agents, these models can also produce biological inco¬ 
herent behaviour of the simulated individuals: the influ¬ 
ence of some congeners can be omitted/overestimated or 
the resulting vector can point towards a direction where 
no stimulus is present (Fig. [T). In addition, such pro¬ 
cess have difficulties in reproducing experimentally ob¬ 
served choices between two concurrent stimuli. For ex¬ 
ample, zebrafish larvae randomly chose to orient towards 
one of two equidistant source of light and do not follow 
the bisector [47]. Here, we present a novel algorithm 
to account for information processing by the individual. 
Rather than summing influential vectors, we propose to 
sum probability distribution functions to orient towards 
the different stimuli. As shown in the model section, such 
mixture distribution can spontaneously produces transi¬ 
tions from choice to consensus and better describes dif¬ 
ferent biological observations. 


D. A multilevel approach 

It would be interesting to develop a multilevel mod¬ 
eling approach that takes into account models for per¬ 
ception and information processing by the individuals. 
Beyond existing SPP and RW models, it requires extend¬ 
ing the description of the agents by biologically relevant 
properties. We propose a new hybrid class of models, at 
the crossroads of both RW and SPP models, that include 
a probabilistic component (inspired by RW) in the be¬ 


FIG. 1: Potential biological incoherence or problematic situ¬ 
ations in the different models assuming a summation of inter¬ 
action vectors between the blue focal fish (FF) at the center 
and green perceived fish (non-perceived fish are represented 
in red). In the metric model, the FF perceives no (A) or a 
high number (B) of neighbors. In the topological model, the 
position of the closest neighbors can be conflicting (C) 
or disproportionate regarding the rest of the group (D). Even 
in visual model, the position of the perceived fish can lead 
to a resulting vector oriented between the stimuli (E) or the 
number of fish can be overestimated by their captured angle 
due to proximity (E). 


haviour of the agents that react to their perception field 
(inspired by SPP models). In the model, the agents chose 
a direction to move according to a probability distribu¬ 
tion function (PDF) that is determined by the presence of 
stimuli in their visual perception field, in accordance with 
observed stochastic choices made by zebrafish larvae m- 
This PDF is computed as a mixture distribution of von 
Mises distributions centered on each perceived stimulus. 
Once the probability distribution function has been com¬ 
puted, the direction of the agent is chosen accordingly to 
this PDF. With this model, we simulate motion of sin¬ 
gle individuals or groups of fish evolving in a bounded 
environment that can include other stimuli like spots of 
interest. 

To validate this model, we measured the individual 
(single fish) and collective (group of 10 fish) motion of 
adult zebrafish in different environmental conditions. Ze¬ 
brafish form loosely cohesive groups that do not show 
strong alignment and could be challenging to model with 
a classical approach. In addition, we observed the be¬ 
haviour of zebrafish in a bounded tank with potential 
spots of interest to take into account the interactions 
with the environment. First, we analysed the locomo¬ 
tion of isolated individual evolving in a uniform environ¬ 
ment to determine the intrinsic characteristics of their 
motion (speed, change in orientation and probability of 
presence). We performed similar measurements in het¬ 
erogeneous environment by placing two floating plastic 
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disk acting as attractive spots. Then, we investigated 
the influence of conspecifics on the spatial distribution of 
fish by observing collective motion of group of zebraflsh in 
both homogeneous and heterogeneous environment. We 
measured their probability of presence in the experimen¬ 
tal tank and their inter-individual distances. This exper¬ 
imental data is used for both, setting parameter values 
for individual behaviour, as well as to compare the pre¬ 
dictions obtained from our model. 


II. NEW STOCHASTIC MODEL 

Our aim is to model fish swimming in nearly 2D in 
a bounded environment with potential spots of interest 
that attract them. To do so, we simulate agents that 
update their position vector Xi with a velocity vector 
Vi though a discrete time process in a bounded two- 
dimensional space: 

Xi{t^5t)=Xi{t)^Vi{t)5t ( 1 ) 

Vi{t + 5t) = Vi{t + 6t)&i(t + 6t) ( 2 ) 

with Vi the linear speed of the agent and 0^ its ori¬ 
entation. Since we focus on the decision making process 
of the fish to chose its orientation in a complex environ¬ 
ment (walls of the tanks, spots of interest) with other fish, 
we assume the simple hypothesis that the linear speed Vi 
of the agent is randomly drawn from the instantaneous 
speed distribution measured in our experiments. The 
novelty of this model is the computation of the orienta¬ 
tion 0^. We model the spherical visual perception field 
of fish (Fig. and describe the probability for fish to 
move in all potential directions by a circular probability 
density function extending from —tt to tt. Thus, 0^ is 
not computed as a resulting vector with potential noise 
but is drawn from a circular probability density function 
(PDF) formed by von Mises distributions, an equivalent 
to the gaussian distribution in circular probability. This 
PDF is characterized by a measure of location /i (equiv¬ 
alent to the mean of a Gaussian PDF) and a measure of 
concentration hi (an inversely proportional equivalent to 
the variance of a Gaussian PDF). For a fish that perceives 
no stimuli, the distribution of probability is described by 
/i = 0 while the value of hcq is determined experimentally 
(Eq.§. By doing so, a fish that perceives no stimuli will 
move forwards with deviation inversely proportional to 
kq. Since our goal is to model fish moving in a bounded 
tank, we introduce the interaction of the fish with the 
walls in the computation of this PDF. As soon as a fish 
is situated in a distance shorter than a threshold value 
dw, we assume that the fish starts a wall-following be¬ 
haviour. To simulate this behaviour, the value of fi is 
not equal to 0 but to the direction along the followed 
wall. Since there are two potential direction, the PDF 
is computed as the weighted sum of the PDF associated 


with each direction. Thus, the PDF /o(6^) for a fish to 
move in each potential direction ^ in a bounded tank 
without perceptible stimuli is given by: 

{ 2'k^q{k.q) ^^p\.'^ocos{e)], \id>dw 

— ^exp[K^cos{0 - ifd<du, 

2 ^ 2TTlo{Kyjj 

(3) 

oo f n\2k 

with d the distance with the closest wall, d^ the thresh¬ 
old distance determining the interaction with the walls, 
kq and Kyj the dispertion parameters respectively asso¬ 
ciated with the basic-swimming and wall-following be¬ 
haviours, fiwi and IJ.W 2 fhe two possible directions along 
a considered wall and Iq the modified Bessel function of 
first kind of order zero. The value of d^^ and are 
determined experimentally. 

Equations Q to Q are sufficient rules to simulate a 
fish swimming in an experimental tank. However, since 
we aim at simulating groups of fish moving in a homoge¬ 
neous but also heterogeneous environment, we implement 
the interactions with other congeners and spots of inter¬ 
est. As soon as the fish perceive stimuli in its perception 
field, its PDF is influenced by those stimuli following two 
steps: information gathering and information processing. 

A. Information gathering 

We simulate fish-agents that can perceive and react 
to 3 dimensional stimuli perceived in their visual field. 
Fish are modeled as 3D polygons with 6 vertices swim¬ 
ming on a 2D plane space (Fig. |2}\). Their visual percep¬ 
tion is simulated as a cyclopean vision sensor that has 
a 270 degree spherical field of view extending frontally 
and laterally and an infinite effective radius (Fig. [^). 
Fish perceive potential stimuli by the solid angle that the 
projection of their extremities capture in their spherical 
perception field (Fig. |^). 

To reflect our experimental conditions, we include two 
potential stimuli in our simulations: fish and spots of 
interest. Fish are considered as polygons of length = 
0.035m, width = 0.01m and height = 0.01m whose ex¬ 
tremities form the two boundaries of the fish. Spots of 
interest are considered as floating disks with a radius of 
0.1m floating 0.05m above the plane space. 

B. Information processing 

Once all potential stimuli have been perceived, the 
model computes a PDF for the focal fish to move ac- 
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FIG. 2: Scheme of the simulated fish and their perception abilities. (A) Simulated fish swim on a 2D plane space and are 
characterized by a width (1 cm), length (3.5 cm) and height (1cm). (B) Their spherical perception held is dehned by an angular 
perception zone and an inhnite distance of perception. The position (x, y) of the hsh is updated at each time step by a velocity 
vector computed according to the stimuli perceived. (C) Fish present in their perception held are perceived as the solid angle 
Af. captured by their extremities. Similarly, we computed the solid angle captured by the spots of interest. 


cording to each stimulus (hsh or spots in this study). 
For example, the probability of the focal hsh to orient 
towards a perceived hsh is given by a von Mises distribu¬ 
tion clustered around this hsh: 


= o . exp[KfCos{0 - fifj] (5) 

Z7rlo[Kf) 

with the potential direction of movement of the hsh, 
/i f. the location of the perceived hsh i and Kf a measure 
of concentration. 

The model draws such distribution for each stimulus 
(hsh or spot) in the perception held of the focal hsh. 
Then for each type of stimuli, it performs a weighted sum 
of all distributions proportional to the ratio of the solid 
angle A^. captures by each stimulus to the sum of the 
solid angles captured by all stimulus At^ . For example, 
the PDF computed for the perceived hsh is given by : 


fpio) = T] 

i=l 


^fi 1 

Atj, 27r/o(K/) 


exp[KfCos{e - Hf,)] 


( 7 ) 


^Tf = ^ Af^ (8) 

i=l 


with Atj. to the sum of the solid angles Af. captured by 
all hsh. Thanks to the summation of PDFs rather than 
vectors, the model intrinsically produces transitions be¬ 
tween consensus (i.e. the hsh orients between two stim¬ 
uli) and choice (i.e. the focal hsh orients toward one of 
two stimuli) according to the angle between to stimuli 
(Fig.||). 

We calculate a similar PDF for the spots of interest 
perceived by the focal hsh: 


fs{0) = E 


^5, 


At^ 27rIo{Ks) 


exp[KsCos{0 — Psi)] 


(9) 


n 

AT, = '^Asi ( 10 ) 

i=l 

with psi the location of the center of the perceived 
spots, t<is the dispersion parameters of the perceived 
spots, Ag. the solid angle captured by the spot and At^ 
the sum of the solid angles A^. captured by all spots. 

Once the PDFs have been computed for each type of 
stimuli (hsh - spots), we calculate a weighted sum of the 
PDFs to obtain the global probability distribution func¬ 
tion f{0) of the focal hsh to move towards a given direc¬ 
tion. In this hrst approach, we assumed that the weight 
of the PDF associated with each type of stimuli is a linear 
function of the total solid angle At^ that they capture in 
the perception held of the focal hsh (Fig.|^. This implies 
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FIG. 3: Evolution of the total probability density function 
friO) according to the position of the neighboring individ¬ 
uals. Contrary to classical models, the fish will not always 
move forward as it would be predicted by weighting the at¬ 
traction vectors. When the neighbors are far from each other, 
the focal fish will randomly orient towards one of them. On 
the contrary, if the two neighbors are closer, the sum of the 
PDF associated with to these fish predict that the focal fish 
will favor the direction between them rather than selecting 
one of them. Thus, by summing the PDFs rather than inter¬ 
action vectors, the model reproduce transitions from choice 
to consensus and vice versa according to the angle between 
two perceived stimuli. 


that the fish respond more strongly to largely perceived 
stimuli but it also allows a potential hierarchy in the re¬ 
sponse to different stimuli. Thus, the global PDF f{0) is 
given by: 



^ + 

Fish fp(0)=y ^ ^ exp[K'fcos(0-M,)] 


A 1 

Spots fs(G)=y exp[K-^cos(6>-/i^.)] 

s 

Total 

-- 


f{0) = 


fo{0) + a^ATffF{0) + /3^AT^fs{0) 


= 



Tf + P^At^ 

0 ^ 0 , 

IV 

dw 

aw, 

if d < 

dw 

Po, 

1^ 

IV 

dw 

Pw, 

if d < 

dw 


( 11 ) 


( 12 ) 


(13) 


with ao and Po the parameters weighting the influence 
of respectively the fish and the shelters for a fish distant 
from any wall and aw and pw for a fish following a wall. 
These parameters are fitted experimentally. 

Then, we numerically compute the cumulated distribu¬ 
tion function (CDF) corresponding to this custom PDF 
f{0) by performing a cumulative trapezoidal numerical 
integration of the PDF in the interval [—7r,7r]. Finally, 
the model draws a random direction 0^ in this distribu¬ 
tion by inverse transform sampling. The position of the 
fish is then updated according to this direction and his 
velocity with equations and 


FIG. 4: Weighted sum of the three PDFs calculated for the fo¬ 
cal fish according to the congeners and the shelters perceived 
by the individual to compute its final PDF. The direction 
taken by a fish is drawn randomly from the final computed 
PDF (in green) by inverse transform sampling. 


III. RESULTS 


A. Homogeneous environment 


We measured the positions of ten fish tested individu¬ 
ally swimming alone in our 1.20m x 1.20m experimental 
tank during 1 hour. Based on this tracking, we build 
the trajectories of the fish and computed their speed and 
change in orientation. An example of a 10 minutes trajec¬ 
tory of a fish is given on Fig. Ik- Fig.[^ shows the cumu¬ 
lated distribution of all instantaneous speeds measured in 
a homogeneous environment with a average speed of 0.07 
± 003 ms“^. Similarly, Fig. represents the cumula¬ 
tive distribution of the changes of orientation along the 
trajectories and highlights that fish are mainly moving 
forwards with soft change of orientation (average change 
= -0.03 ± 0.84 rad). The distribution of the positions 
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detected in the tank (Fig. [n^v) shows that the higher 
probability of presence was found along the walls. Thus, 
in a homogeneous environment, fish were mainly swim¬ 
ming along the walls of the tank and avoid the centre of 
it. 
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FIG. 5: Example of experimental (A-B) and simulated (C-D) 
trajectories of a fish swimming alone during 10 minutes in the 
absence (A-C) or presence (B-D) of two floating disks. The 
colour of the trajectory indicates the speed of the individual. 
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We used this experimental data to set the parameters 
of our model in order to simulate the movement of a sin¬ 
gle fish in a homogeneous tank. In the simulations, the 
speed of the agent is drawn from the experimental dis¬ 
tribution of the instantaneous speed of the fish. Speeds 
are drawn independently from each other so that there is 
no correlation between the speed of agent i at time t and 
t + 1. While this differs obviously from the reality, for 
simplicity we do not take into account speed matching 
in this first parametrization of our new model. Experi¬ 
ments with single fish allow us to fit the parameter value 
characterizing the change of direction of fish hq. To do 
so, we measured the change of direction of the fish when 
they were at least at 0.30m from any walls. By doing 
so, we excluded the potential influence of the walls and 
consider only the intrinsic change of direction. Then, we 
compared this experimental distributions to theoretical 
ones and select the best fitting. In this first approach, we 
chose the sum of least squares minimization as a rapid 
and low computational costly fitting method. The best 
fitting of this experimental distribution was obtained by 
/^o = 6.3, as shown by Fig. 

For experiments involving a single fish, only the inter¬ 
action with the walls of the tank is present. Therefore, 
the other relevant parameters are the distance of interac¬ 
tion with the walls dyj and the measure of concentration 
Kw of the PDFs associated with wall following. To esti¬ 


FIG. 6: Experimental individual behaviour of single fish from 
AB strain in the experimental tank without stimuli. (A) The 
distribution of the speed shows an average speed of 0.07 ± 
0.03 ms“^. (B) The distribution of the change in orientation 
highlights that fish are mainly swimming forward with low 
deviation. Results are cumulated for 10 replicates of one hour. 

mate these parameter values, we performed simulations 
with different couples of value k^j) and compared 
the experimental distributions of change in orientation 
and of probability of presence with those generated by 
the simulations. Fitting of these parameters showed that 
dw = 0.05 and = 20 are the best values to reproduce 
our experimental data (Fig. [^. A 10 minutes trajectory 
of a simulated fish with these parameter values is shown 
on Fig[^. 

Then, we performed experiments with 10 groups of 
10 zebrafish. As observed for single fish experiments, 
fish were mostly detected along the walls of the tank 
(Fig. prjG). Since we did not track the fish individually 
when swimming in group, we could not build the indi¬ 
vidual trajectories of the fish. However, we measured 
the distance between all pair of fish at each time step as 
a measure of group cohesion. The distribution of these 
inter individual distances shows an average of 0.394m ± 
0.38m with the mode of the distribution between 0.1 and 
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FIG. 7: Experimental distribution and fitting by sum of least 
squares minimisation of direction changes of single zebrafish 
swimming in a homogeneous tank. (A) Experimental dis¬ 
tribution (blue) and fitting (red) of changes of direction ob¬ 
tained for = 6.3. The data were calculated for trajectories 
measured at least at 0.3m from any walls to minimize their 
influence on the fish orientation. (B) Evolution of the sum of 
the mean squares between the experimental distribution and 
a theoretical distribution of 100,000 random draws in a von 
Mises distribution with parameters /x = 0, kq. 


0.2m (Fig. 1^). 

To simulate experiments with 10 fish, we introduce 
three parameters describing the interaction with the fish: 
the measure of concentration hip associated with the PDF 
computed for each fish and the parameters and aw 
that weight the influence of other individuals on a fish 
that is far from (ao) or close (aw) to a wall. The value 
of is assumed to be similar than k.w and is equal to 
20. By doing so, we consider that a fish orients towards 
a given target /i with a high accuracy. To determine the 
value of the weights we perform simulations with dif¬ 
ferent values of (o^o, aw) and compare the distribution of 
interindividual distances and probability of presence with 
those obtained for the experiments. The best values to 
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FIG. 8: Fitting by sum of least squares minimisation of the 
parameters dw and kw determining the interaction of a single 
fish with the walls of the tank. The parameter dw corresponds 
to the threshold distance of interaction with the walls and 
Kw is inversely proportional to the width of the von Mises 
distribution associated with wall following. To each couple 
(dw, i^w) corresponds a couple of squares whose colours in¬ 
dicate the value of the sum of least squares obtained for the 
comparison of the changes of direction (left square) and the 
probability of presence (right square) with the experimental 
data shown in Fig. and Fig. [nj\. Lower values of sum of 
least squares corresponding to better fits are in the purple and 
blue colour scale. (A) The exploration of a first set of param¬ 
eter values indicated that the best couple {dw and kw) are 
found for values {dw = 0.05 and 15 < kw < 30), highlighted 
in the red rectangle. (B) The refinement of the parameters 
exploration in the ranges determined in (A) showed that the 
best fitting is obtained by dw = 0.05 and kw = 20. 


reproduce both distributions are ao = 55 and aw = 20 
(Fig.[^. Thus, fish that follow a wall are less influenced 
by other congeners than fish situated in the center of 
the tank. With these parameter values, our model was 
able to reproduce the probability of presence displayed 
by groups of 10 zebrafish as shown in Fig. pTp. Con¬ 
cerning the distribution of interindividual distances, the 
model reproduces the decreasing distribution observed in 
the experiments except for the mode of the distribution 



































that is between Om and O.lm (Fig. [^). 




FIG. 9: Cumulated interindividual distances measured be¬ 
tween all pairs of fish. Results are obtained for 10 experiments 
of one hour with 10 groups of 10 fish (A) and 10 simulations 
with 10 agents (B). 
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FIG. 10: Fitting by sum of least squares minimisation of the 
parameters ao and aw determining the interaction of a fo¬ 
cal fish that far away or close to a wall with other individu¬ 
als. To each couple (ao, olw) corresponds a couple of rect¬ 
angles whose colours indicate the value of the sum of least 
squares obtained for the comparison of the interindividual dis¬ 
tances (left rectangle) and the probability of presence (right 
rectangle) with the experimental data shown in Fig. and 
Fig.[njB. Lower values of sum of least squares corresponding 
to better fits are in the purple and blue colour scale. (A) The 
exploration of a first set of parameter values indicated that 
the best couples (ao, aw) are found for values (ao < 100, 
10 < aw < 100), highlighted in the red rectangle. (B) The 
refinement of the parameters exploration in the ranges de¬ 
termined in (A) showed that the best fitting is obtained by 
ao = 55 and aw = 20. 


We added two spots of interest in the experimental 
tank placed at 25v^ cm from two opposite corners along 
the diagonal of the tank. These spots consisted of blue 
plastic disks (20cm of 0) floating at the water surface and 
hung by nylon threads. An example of path followed by 
a fish during 10 minutes is shown on Fig. [^. We calcu¬ 
lated similar parameters from the individual positions of 
fish (10 X 1 fish) moving alone in the tank with two spots. 
In the presence of two spots, the fish are mainly detected 
along the walls and under the spots as shown by their 
probability of presence (Fig.[T^). Experiments with sin¬ 
gle fish also showed that fish decreased their speed under 
these spots. Indeed, the separation of the speed distri¬ 
bution measured Outside or Inside the spots shows that 
the average speed was three times slower when fish were 
located under a floating disk (Fig.jl^). On the contrary, 
the presence of the floating disks did not affect the dis¬ 
tribution of the changes in orientation that were similar 
Outside and Inside (Fig. [T^). 


For experiments involving a single fish, only the inter¬ 
action with the walls and the shelters of the tank are 
present. Since d^ = 0.05 and = 20 were fitted by our 
experiments in homogeneous environment, we explore the 
parameter values of /3q and flw (ihe ponderation of the 
influence of the spots) and assume that the measure of 
concentration associated with the spots ns = i^w = 20. 
We perform simulations for different values of (/3o,/dvp) 
and compare their results with the experimental proba¬ 
bility of presence and distribution of change of orienta¬ 
tion. The best value to reproduce our experimental data 
are l3o = 0.15 and Pw = 0.01 as shown by Fig. As 
observed for the influence of the congeners, these values 
indicate that fish following a wall are less influenced by 
the spots than fish swimming in the centre of the tank. 
With these parameter values, the model reproduce the 
spatial distribution of the fish along the wall and under 
the floating disks (Fig. [T^). 

Finally, we observed 10 groups of 10 zebrafish swim- 
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FIG. 11: Probability of presence for experimental (A, B) and 
simulated (C, D) data. Results are obtained for 10 replicates 
of Ih. (A) Experimental individual behaviour of single fish 
from AB strain in the empty experimental tank. The prob¬ 
ability of presence showed that fish were mainly swimming 
along the walls of the tank. (B) Experimental results ob¬ 
tained for groups of 10 AB strain zebrafish in a homogeneous 
environment. Similarly to single fish, groups of 10 zebrafish 
were mainly swimming along the walls of the tank. (C) Prob¬ 
ability of presence of a simulated single fish in a homogeneous 
environment during one hour (results for 10 simulations). (D) 
Probability of presence of a simulated group of 10 fish in a 
homogeneous environment during one hour (results for 10 sim¬ 
ulations) . 


ming in the presence of two spots of interest during one 
hour. In this case, the fish were also observed mainly 
under the spots and along the walls but show a prefer¬ 
ence for the spots (Fig. [T^). The measure of the inter¬ 
individual distances shows that the presence of spots does 
not have a strong influence on the distance between the 
fish (Fig.[T^). Indeed, the average distance between the 
individuals is 0.439m ±0.331m, which is 6mm less than 
observed in the absence of floating disks. 

We simulated these experiments with groups of 10 ze¬ 
brafish with our model by integrating the interactions 
of the agents with both the other fish and the spots of 
interest. In the previous experiments, the parameter rul¬ 
ing the interactions with the fish (o^o, o^w) and the spots 
of interest (/3o, Pw) were fitted independently but here, 
both stimuli are simultaneously present in the percep¬ 
tion field of the fish. Therefore, we investigate the rel¬ 
ative importance of both stimuli by weighting the influ¬ 
ence of the other fish or the spots of interest. To do so, 
we multiply the previously fitted values of (o^o, olw) and 
(/3o, Pw) by different weighting factors wp = 
and ws = example, values wp = \ and 

ws = \ imply that the two stimuli are simply added 
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FIG. 12: Experimental individual behaviour of single fish 
from AB strain in the experimental tank with two spots of 
interest. (A) The distribution of the speed shows an aver¬ 
age speed of 0.09 ± 0.06 ms“^ outside the spots while fish 
were swimming with an average speed of 0.03 ± 0.02 ms“^ 
under the spots. (B) The distribution of the change in orien¬ 
tation highlights that fish are mainly swimming forward with 
low deviation both outside and inside the spots. Results are 
cumulated for 10 replicates of one hour. 


while values wp = ^ and ws = | imply that the influ¬ 
ence of the fish is divided by 2 and the influence of the 
spots by 5. Thus, we perform simulations for each couple 
of {wp^ws) and compare the distributions of the proba¬ 
bility of presence and the interindividual distances with 
those measured experimentally. The best fit is given by 
wp = 2 and ws = 9 (Fig.p!^ which implies that the influ¬ 
ences of the fish and the spots have to be decreased by a 
factor 2 and 9 respectively. Therefore, when both stimuli 
(congeners and spots) are perceived by the fish, the best 
values of a* and to reproduce the experimental data 
are = 27.5, aw = 10, /3o = 0.016 and Pw = 0.0011. 
While these values give a correct fit of the experimen¬ 
tal probability of presence (Fig. [T§d) , the model do not 
perfectly reproduce the distribution of the interindivid¬ 
ual distances (Fig.[T^). Similarly than for the groups of 
fish in a homogeneous environment, the mode interindi- 
vidual distances measured in the simulation is between 0 
and 0.1m while the mode of the experimental distances 
is between 0.1 and 0.2m. 
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FIG. 13: Probability of presence for experimental (A, B) and 
simulated (C, D) data. Results are obtained for 10 replicates 
of Ih. (A) Experimental individual behaviour of a single fish 
from AB strain in a tank with two floating discs. The prob¬ 
ability of presence showed that fish were mainly swimming 
along the walls or under the discs. (B) Experimental results 
obtained for groups of 10 AB strain zebrafish in a heteroge¬ 
neous environment. The probability of presence showed that 
fish were mainly present under the floating disks. Results 
are obtained from 10 replicates of one hour. (C) Probability 
of presence of a simulated single fish in a homogeneous en¬ 
vironment during one hour (results for 10 simulations). (D) 
Probability of presence of a simulated group of 10 fish in a 
homogeneous environment during one hour (results for 10 sim¬ 
ulations) . 


IV. DISCUSSION 



EIG. 14: Pitting by sum of least squares minimisation of the 
parameters /3o and determining the interaction of a focal 
fish that far away or close to a wall with the spots of interest. 
To each couple (/3o, ^w) corresponds a couple of rectangles 
whose colours indicate the value of the sum of least squares 
obtained for the changes of orientation (left square) and the 
probability of presence (right square) with the experimental 
data shown in Eig.[T^ and Eig. [T^. Lower values of sum of 
least squares corresponding to better fits are in the purple and 
blue colour scale. (A) The exploration of a first set of parame¬ 
ter values indicated that the best couples (/3o, fiw) are found 
for values (0.05 < /3o < 0.5, 0 < < 0.25), highlighted 

in the red rectangle. (B) The refinement of the parameters 
exploration in the ranges determined in (A) showed that the 
best fitting is obtained by /3o = 0.15 and fiw — 0.01. 


A. Experimental results 


In this study, we observed individual and collective be¬ 
haviour of zebrafish in homogeneous and heterogeneous 
environments. Firstly, we observed single D. rerio swim¬ 
ming in our experimental tank to describe the swimming 
behaviour of zebrafish. In homogeneous environment, the 
fish were mainly following the wall of the tank and avoid 
the center of it. This observation was also reported in 
studies performed with zebrafish and robotic-fish [48]. 
We characterized the motion patterns of single individ¬ 
uals and found out an average speed of 0.07 ±0.03 m/s 
and an average change of orientation of —0.03 ±0.84 rad, 
in coherence with previous studies performed on single 
zebrafish individuals [31 SOI HU- Our experiments with 
10 zebrafish showed that the spatial repartition of the 
group did not differ significantly from the individual one. 
Fish in groups were also mainly detected along the walls 
of the tank. 

The presence of floating disks influence the spatial 


distribution of the fish in the experimental tank. The 
shaded area seemed as attractive as the walls since fish 
showed a similar probability of presence for both stimuli. 
These disks had no influence on the change of direction 
of the fish but resulted in a spatial differentiation of the 
instantaneous speed of the individuals. Indeed, the fish 
swam with an average speed of 0.09 ± 0.06 m/s outside 
the shelters but with an average speed of 0.03 ± 0.02 
m/s under them. Such reduction of speed can indicate 
that shaded areas are potentially considered as tempo¬ 
rary resting sites by the fish. The observation of groups 
of 10 fish showed similar results with preference for both 
the walls and the shaded area. 

These results show that zebrafish are avoiding free wa¬ 
ter and prefer to swim near potential shelters (floating 
objects or bank). Although fish were attracted by the 
floating disk, they did not seem to avoid the light since 
their presence under the disk and along the wall (that are 
exposed to light) are similar. These observations are in 
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FIG. 15: Fitting by sum of least squares minimisation of the 
parameters weighting the influence of the fish wp and the 
spots ws in simulation with 10 zebrafish swimming with two 
spots of interest. To each couple {wf, ws) corresponds a cou¬ 
ple of rectangles whose colours indicate the value of the sum 
of least squares obtained for the comparison of the interindi¬ 
vidual distances (left square) and the probability of pres enc e 
(right square) with the experimental data shown in Fig. 16 \ 
and Fig.p^^. Lower values of sum of least squares correspond¬ 
ing to better fits are in the purple and blue colour scale. The 
exploration of a set of parameter values indicated that the 
best couples {wf, ws) is icf = ws — 


accordance with the ecology of the species that shows a 
diurnal activity in the Nepalese and Indian shallow water 
and rice paddies m- 




FIG. 16: Gumulated interindividual distances measured be¬ 
tween all pairs of fish swimming in the presence of two spots 
of interest. Results are obtained for 10 experiments of one 
hour with 10 groups of 10 fish (A) and 10 simulations with 
10 agents (B). 


B. A new hybrid modeling approach 

Based on our observations and on the literature, we de¬ 
veloped a model describing individual and collective mo¬ 
tion. Currently, collective motion receives attention from 
researchers of numerous fields such as biology, physics, 
computer science and robotics. Current modeling tends 
to give up traditional methods (metric or topologic) ac¬ 
counting for the perception of influential neighbors by 
the focal fish for vision-based approach j43l445] . Here we 
present a 3D perception system accounting for zebrafish 
vision in which the objects are described by the solid an¬ 
gle that they capture in their perception field. By doing 
so, we are closer to a realistic description of the sensory 
system of the fish. In this first step, we assumed the sim¬ 
ple hypothesis that object were homogeneously perceived 
in this 3D sensory field. However, the recent characteri¬ 
zation of the visual system of zebrafish highlighted cen¬ 
ters of acute vision (i.e. areae) in the fronto-dorsal region 
[46] . Therefore, future models should take into account 
such heterogeneity in the perception field. In addition, 
the position of the areae differs from one species to the 


other and is expected to influence the structure of the fish 
school [46] . The extension of our model to 3 dimensional 
movements would allow to test such hypothesis. 

In addition to this sensory system, we proposed a new 
mechanism to determine the direction of an agent accord¬ 
ing to its visual field. Rather than computing a resulting 
vectorial force that is applied to the fish, we redefined 
the modeling approach by describing the decision-making 
process of the fish to make an intentional movement ac¬ 
cording to its perception field. In this first step, we rep¬ 
resented the choice of the individual by a stochastic pro¬ 
cess. The fish can potentially move in any direction but 
it will favor directions associated with a perceived stim¬ 
uli (congeners for example). This stochastic model can 
reproduce collective behaviour exhibited by species that 
show cohesive behaviour without presenting higher order 
or only occasionally without introducing a high level of 
noise accounting for disturbances. This is made possible 
by the representation of the perception field of the fish 
and its translation into a probability distribution func¬ 
tion. By doing so, we are closer to an effective description 
of the individual decision-making process during motion 
[46] . Indeed, with such approach of PDFs’ summation. 
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we can include potentially all kinds of stimuli (congeners, 
environment, food...) that are perceived by the fish and 
account for a choice or consensus between potentially an¬ 
tagonist stimuli. Such choice between two concurrent 
stimuli has been evidenced in zebrafish larvae that orient 
towards one of two light source rather than swimming 
towards their bisector 07]. The authors showed that two 
retinal pathways controlling turn movements and rapid 
forward swimming are responsible for phototaxis. Such 
advances on the understanding of information process¬ 
ing by neuronal pathways should be integrated in further 
multi-scale models to make the link between fish move¬ 
ment, perception of information and its processing based 
on biological knowledge. 

Moreover, most of the models consider motion in un¬ 
bounded homogenous space (torus geometry). While this 
hypothesis is reasonable to study the collective behaviour 
of animal living in pelagic water, the interactions with 
can not be neglected for species living in small streams. 
We extend our model to take into account bounded and 
heterogenous space getting closer to natural conditions of 
fish like D. rerio. The model correctly reproduced the be¬ 
haviour of single individual swimming in a bounded tank. 
The fish mainly follows the wall but sometimes swim in 
the center of the tank. The probability of presence was 
also correctly reproduce but the distribution of interindi¬ 
vidual distances was biased towards short distances. This 
could originate from the absence of preferred distance or 
avoidance distance between the agents. Indeed, while the 
agents of the model can overlap, the fish have to respect 
a distance between them. Such distance was not intro¬ 
duce in this first parametrization of the model to limit 
the number of parameters. In a heterogeneous environ¬ 
ment, the model was also able to reproduce the individual 
trajectories of single fish that transits between the walls 
and the spots of interest. As in the homogeneous envi¬ 
ronment, the probability of presence of groups of fish is 
also fitted, but again, the distribution of the interindivid¬ 
ual distances is biased towards short distances. There¬ 
fore, these results indicate that the approach developed 
in this study can successfully describe individual and col¬ 
lective motions of fish but that further version could in¬ 
clude additional biological variable and behaviour that 
are species dependent. For example, here the probability 
distribution function for a simulated fish are centered on 
the position of the other fish that it perceives but the 
PDF could be centered on the direction of the congeners 
to simulate an alignement behaviour. Rather than de¬ 
veloping a model that exactly described the behaviour 
of zebrafish, our goal was to introduce a new decision¬ 
making algorithm for motion in complex environment. 

In this study, we applied our decision making algorithm 
to an individual-based model (IBM). Recently, a kinetic 
model (KM) was proposed to account for the individual 
movement of isolated zebrafish 03110 ]. This model in¬ 
spired by [37l [38] took into account a dynamic speed reg¬ 
ulation characterizing D. rerio motion. In our model, we 
did not implement a function accounting for speed mod¬ 


ulation in order to focus on the decision-making mecha¬ 
nism. The instantaneous speed was drawn from the ex¬ 
perimental distribution with a time step corresponding 
to the tail-beat period of zebrafish. Then, future work 
should investigate the impact of the proposed stochastic 
mechanism for orientation on such KMs. Moreover, the 
development of a KM version of our IBM could be an 
intermediate step towards a continuum model (CM) de¬ 
scription of the group of agents. This could then be use to 
perform large-scale analysis and prediction of the collec¬ 
tive behaviour displayed by very large population. Such 
multiscale modeling approach would allow us to identify 
the properties of the group that are preserved at all scales 
of analysis or on the contrary that are specific to a partic¬ 
ular level of observation [51]. We performed simulations 
with a small number of agents to mimic our experimen¬ 
tal conditions in this study that are close to the natural 
size of zebrafish group [50] . While our IBM could model 
larger groups of tens of individuals, simulations involving 
thousands of agents would require a longer computational 
time than KM and CM. 


In parallel to the understanding of information pro¬ 
cessing by individuals, this approach is also a new step 
towards bio-inspired algorithms that can be implemented 
in robotic agents. Indeed, it is a major scientific challenge 
to build artificial systems composed of robots that can 
perceive, communicate to, interact with other agents (bi¬ 
ological or robotic) and adapt to their environment [52]. 
To do so, we need to develop artificial agents that com¬ 
municate through appropriate channels corresponding to 
specific animal traits but also that correctly perceive and 
interpret signals emitted by the animals [SllM]. This 
was firstly achieved in 2007 by building bio-inspired ar¬ 
tificial cockroaches that where able to sense the presence 
of congeners and to adapt their behaviour following a 
bio-inspired algorithm [55] . 


In fish, an increasing number of studies aim at devel¬ 
oping such robotic agents to interact with group of fish 
[56ll59] . Current experiments generally involve one robot 
that follows a predetermined trajectory or that moves ac¬ 
cording to the position of fish detected through the inter¬ 
mediary of a camera that films the entire setup. While 
this methodology is in the line of the classical ethological 
experimentation to investigate behavioural responses of 
animals, the development of fully autonomous integrated 
lures in fish schools (or other species) requires the devel¬ 
opment of perception abilities for the robotic agents as 
well as adapted behavioural algorithms. 


In this perspective, the development of perception- 
based models is a necessary step towards intelligent ar¬ 
tificial systems capable of closing the loop of interaction 
between animals and robots. 
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V. METHODS 

A. Animals and housing 

We acquired 150 adult wild-type zebrafish {Danio re- 
rio AB strain from Institut Curie (Paris). Fish were 18 
months old at the time of the experiments. We kept fish 
under laboratory conditions, 27°C, 500//S salinity with 
a 10:14 day might light cycle. The fish were reared in 
55 litres tanks and were fed two times per day (Special 
Diets Services SDS-400 Scientific Fish Food). Water pH 
was maintained at 7.5 and Nitrites (NO^“) are below 0.3 
mg/1. 

B. Experimental setup 

We recorded the behaviour of zebrafish in a 120 x 120 
X 30 cm experimental tank made of glass with internal 
walls covered with white adhesive. The water depth was 
kept at 10cm during the experiments. A Logitech®HD 
Pro Webcam C920 was mounted 160 cm above the water 
surface to record experiment at a resolution 1920 x 1080 
and at 15 FPS. The camera was connected to a worksta¬ 
tion Dell®Precision T5600 dedicated to the recording of 
the videos and their analysis. One halogen lamp (450W) 
was placed at each corner of the tank and oriented to¬ 
wards the wall to provide indirect lightning of the tank. 
The whole set-up was confined behind white sheets to 
isolate experiments and homogenize luminosity. 

C. Experimental procedure 

We recorded the behaviour of zebrafish swimming in 
our tank during one hour. We tested two numbers of in¬ 
dividuals (single fish or groups of 10 fish) in two environ¬ 
mental conditions (homogeneous or heterogeneous). The 
heterogeneous environment was created by adding two 
floating disks made of blue coloured Plexiglas (200 mm 
diameter and 3 mm thick) suspended by nylon threads. 
The two spots were spaced from 70 cm and located on a 
diagonal of the square. Before the trials, fish were placed 
with a hand net in a cylindrical arena (20 cm diameter) 
made of Plexiglas placed in the centre of our experimen¬ 
tal tank. Following a five minutes acclimatization period, 
this arena was removed and fish were able to freely swim 
in the experimental tank. We performed 10 replicates for 
each combinaison of parameters (number of fish x envi¬ 
ronmental condition). 

D. Data analysis 

The recorded videos were analyzed offline using a cus¬ 
tom Matlab script developed to detect the position of 
the fish. This script performed a background subtrac¬ 
tion on each frame and transformed it in a binary image 


according to a pixel intensity threshold given by the user. 
The software then identified the blobs on this image and 
kept only the blobs that were formed by more than 20 
and less than 200 pixels (those values were obtained by 
manually identifying the fish on multiple frames). Since 
this method did not allow a perfect detection of all the 
individuals, we developed a second script that was run 
after the first one and that plotted the frame where a 
fish (or more) was undetected for the user to manually 
identify the missing individual(s). While this analysis 
tool is time-costly, it allowed us to identify the fish that 
were partially hidden during a collision/superposition 
with another fish or the fish that were situated under 
the floating disks since our program was not able to de¬ 
tect them due to a lack of intensity of the pixel after the 
background subtraction. The positions P(x, y) of the 
fish were recorded at each time step t = 1/lbs during 
the experiment with single fish in homogeneous environ¬ 
ment and t = l5 for all other experimental conditions. 
This allowed us to build the trajectory of each individ¬ 
ual for experiments involving a single fish and to com¬ 
pute the speed of the individuals as well as their change 
in orientation between successive positions. The instan¬ 
taneous speed Vt was calculated on three positions and 
thus computed as the distance between P(x, t — 1) and 
P(x, t -h 1) divided by 2 time steps while change in ori¬ 
entation were computed as the angle between two succes¬ 
sive speed vectors). The distributions of speed were com¬ 
puted only for parts of the trajectory during which the 
fish were not in freezing behaviour (i.e. immobile). This 
corresponds to a spontaneous speed higher than 1mm 
per second. Since our tracking system did not resolve 
collision with accuracy, we did not calculate individual 
measures for data obtained with groups of fish but char¬ 
acterized the aggregation level of the group. 


E. Implementation and numerical simulations 

The model was implemented in a Matlab script. The 
simulations were run during 10800 time steps, each time 
step representing an increment of 1 /3 second to simulate 
a total time of 1 hour, similarly to our experiments. This 
time step was chosen according to the tail beat frequency 
of the zebrafish of 2.5 Hz. By doing so, we assume that 
zebrafish can potentially change their orientation at each 
tail beat. The position of the agents on the 2D plane 
space is described by the position of their head (x, y, 0). 
The positions of the other vertices are computed accord¬ 
ing to the position of the head and the direction of the 
fish. To compute the solid angle of the perceived stim¬ 
uli in the perception held of the focal fish, we calculated 
the area of the spherical polygon delimited by the pro¬ 
jection of the stimuli’s vertices on a unit sphere centered 
on the focal fish. To do so, we divided the polygon in 
two spherical triangle and computed their spherical ex¬ 
cess using L’Huilier’s theorem : 
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= Y tan(is) tan[i(s - a)] tan[i(s - b)] tan[i(s - c)] 

(14) 

with a, 6, and c the length of the arcs between the vertices 
expressed in spherical coordinates (p, 0^ (j)) and computed 
by: 

a = arccos(sin(0i) sin((/) 2 ) + cos(^i — O 2 ) cos(^i) cos(^ 2 ) (15) 
and s the semi parameter given by : 

a c 

s = ( 16 ) 
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